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ABSTRACT 

We have studied the dependence of topology of large scale structure on tracer, gravitational evo¬ 
lution, redshift space distortion, and cosmology. A series of large N-body simulations of the ACDM 
and SCDM models that have evolved 1.1 or 8.6 billion particles, are used in the study. Evolution of 
the genus statistic, used as a topology measure, from redshift 8 to 0 is accurately calculated over a 
wide range of smoothing scales using the simulations. The tracers of large scale structure considered 
are the CDM matter, biased peaks in the initial density field, dark halos, and ‘galaxies’ populating 
the dark halos in accordance with a Halo Occupation Distribution (HOD) model. 

We have found that the effects of biasing, gravitational evolution, and initial conditions on topology 
of large scale structure are all comparable. The redshift space distortion effects are relatively small 
down to about 5 /i _1 Mpc for all tracers except for the high threshold part of the genus curve. The 
gravitational effects are found to be well-modeled by the analytic perturbation theory when the CDM 
matter distribution is considered. But the direction of gravitational evolution of topology can be 
even reversed for different tracers. For example, the shift parameter of the genus curve evolves in 
opposite directions for matter and HOD ‘galaxies’ at large scales. At small scales there are interesting 
deviations of the genus curve of dark halos and ‘galaxies’ from that of matter in our initially Gaussian 
simulations. The deviations should be understood as due to combined effects of gravitational evolution 
and biasing. This fact gives us an important opportunity: topology of large scale structure can be 
used as a strong constraint on galaxy formation mechanisms. At scales larger than 20 /i -1 Mpc all 
the above effects gradually decrease. With a good knowledge of the effects of non-linear gravitational 
evolution and galaxy biasing on topology one can also constrain the Gaussian random phase initial 
conditions hypothesis to high accuracy. 

Subject headings: 


1. INTRODUCTION 

Topology analysis has long been applied to galaxy red¬ 
shift survey data to test one of the major predictions of 
simple inflationary scenarios: that the primordial density 
fluctuations are Gaussian random phase (Gott, Melott, 
& Dickinson 1986; Gott et al. 1989; Vogeley et al. 1994). 
To recover the primordial density fluctuations it is nec¬ 
essary to smooth the observed sample over large enough 
scales to reach the linear regime. Since redshift surveys 
have not been very deep in the past, one was left with a 
small dynamic range for topology study after big smooth¬ 
ings, and it has been hard to draw firm conclusions on 
the nature of the initial density field. As ambitious new 
redshift surveys like the 2dF and the Sloan Digital Sky 
Surveys)http://sdss.org) are being completed, it is now 
hoped that the topology analysis can become one of the 
precision measures for cosmology. These surveys are not 
just large in angle and deep in depth, but are also dense. 
This means that the new survey samples can be used 
to study the primordial density field at large scales, and 
also to explore small-scale phenomena like the formation 
of galaxies. 

Since the main purpose of studying topology has been 
to discover the Gaussianity of the linear density held, 
there has been relatively little work on topology at small 


non-linear scales. One work on small scale topology is 
Springel et al. (1998) who applied the analysis to their 
N-body simulations and the 1.2-Jy IRAS galaxy sample. 
However, without quantitatively studying the complete 
scale- and cosmology dependences of topology, as mea¬ 
sured by the genus, they have incorrectly concluded that 
the shape of the genus curve is that of a random phase 
held far into the non-linear regime, and that the ampli¬ 
tude drop of the genus curve due to the phase correlation 
is the only sensitive measure of non-linearity. In this pa¬ 
per we will show that topology analysis becomes useful 
when a sample is studied over a wide range of smooth¬ 
ing scale and when quantitative measures of the devia¬ 
tion of the genus curve from the Gaussian one are used. 
In particular, we will present our finding that topology 
of underdense regions at small scales is sensitive to the 
structure formation mechanism. This fact opens the pos¬ 
sibility that galaxy formation can be better understood 
by studying the small scale topology of galaxy distribu¬ 
tion. 

Previous analyses of observational samples for topol¬ 
ogy study have generally used the dark matter distri¬ 
bution for comparison. The similarity of the observed 
genus curves to those obtained from simulated matter 
distribution has been investigated for various cosmogo¬ 
nies. This may be fine when the uncertainties in genus 
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measured from observations are much larger than the ef¬ 
fects of biasing and redshift space distortion. With the 
unprecedented accuracy provided by observational data 
like the SDSS sample, however, these effects have to be 
accurately accounted for. In this paper we will study the 
gravitational evolution of the genus of peaks in the initial 
density field and dark halos as well as dark matter to see 
the dependence of topology on these tracers. Redshift 
space distortion effects are also studied. They all have 
non-negligible, important effects on the genus curve. 

Our work is numerical since the main purposes of our 
work is to explore the strong non-linear effects of gravi¬ 
tational evolution through topology analysis where cur¬ 
rent analytic approaches are inadequate. Study of the 
tracer dependence of, and biasing effects on, topology 
also requires numerical simulations. To check for the im¬ 
portant question of whether the topology is consistent 
with Gaussian random phase initial conditions to high 
accuracy in the new large observational data sets given 
that the effects of non-linear evolution and biasing are 
detectable in such surveys, the best method would seem 
to be direct comparison with large cosmological simu¬ 
lations with Gaussian random phase initial conditions 
where non-linear evolution and biasing are properly mod¬ 
eled. If agreement is obtained, this supports the Gaus¬ 
sian random phase initial conditions hypothesis to high 
accuracy. 

Introduction to the genus and its related statistics is 
given in sections 2 and 3. Our N-body simulations used 
in this study are described in section 4. Evolution of 
the genus curve is studied in section 5 for matter and 
‘galaxy’ distributions in the ACDM and SCDM cosmogo¬ 
nies. Conclusions follow in section 6. 

2. GENUS 

To measure topology of a tracer distribution we cal¬ 
culate the genus of iso-density surfaces as a function of 
density threshold level. The genus is defined as 

G = Number of holes in contour surfaces 

— Number of isolated regions (1) 

in the isodensity surfaces at a given threshold level. In 
this definition, a sphere has genus of —1, and a torus has 
genus 0 (Gott, Melott & Dickinson 1986). The Gauss- 
Bonnet theorem connects the global topology with an 
integral of local curvature of the surface S, i.e. 

G = — — [ KdA, (2) 

4tt Js 

where n is the local Gaussian curvature. The genus 
is equal to —1/2 times the Euler characteristic, or the 
fourth Minkowski Functional in 3-dimensions (Mecke, 
Buchert, & Wagner 1994; Schmalzing & Buchert 1997). 

For a Gaussian random field an analytic formula exists 
for the genus-threshold density relation (Hamilton, Gott, 
& Weinberg 1986). The genus per unit volume is 

g{v) = A{\ - ^e-" 2 / 2 , (3) 

where v = (p — p)/a is the threshold density in unit of 
standard deviations a = ((p — p ) 2 ) 1 / 2 from the mean. 
The amplitude is 

( 4 ) 


where the second moment 

(k 2 ) = j P(k)W(k)k 2 d 3 k/ J P{k)W{k)d 3 k (5) 

depends only on the shape of the power spectrum and 
the smoothing kernel W. Instead of using v to iden¬ 
tify isodensity contours we use the volume fraction label 
Vf at which the corresponding isodensity contours en¬ 
close the fraction of the sample volume equal to that 
at the density threshold Vf<j in the case of a Gaussian 
field. In this way the genus curve becomes independent 
of the one-point probability distribution function, which 
can be better studied by different statistics. Use of Vf 
also makes the genus curves less sensitive to details of 
bias (Park & Gott 1991), which is beneficial in studying 
the primordial field. 

3. GENUS-RELATED STATISTICS 

3.1. Shift Parameter 

Any departure of the genus curve from the relation 
given by equation (3) indicates non-Gaussianity of the 
field. To quantify the deviation we define a set of genus- 
related statistics. The first statistic Av measures the 
shift of the middle part of the genus curve. We compute 
Av defined by (Park et al. 1992) 

Av = J vG ohs (v)dv/ J Gfn(v)dv, (6) 

where Gobs(^) is the measured genus and Gr t(v) is the 
random-phase curve best fit to the measured genus data 
from v = —1 to +1. The best fitting amplitude A 0 b s of 
the observed genus curve is obtained from Gg t (0). 

3.2. Amplitude Drop 

The second statistic is the amplitude drop of the genus 
curve. The density field with phase correlation due to 
gravitational evolution typically has a fewer number of 
structures compared to that with the same power spec¬ 
trum but with random phases. A derived statistic is 

Ra = Aobs/Aps, (7) 

where A P g is the amplitude expected for a Gaussian field 
which has the power spectrum (PS) of the evolved parti¬ 
cle/‘galaxy’ distribution. Existence of phase correlation 
reduces Ra below 1 (Vogeley et al. 1994; Canavezes et 
al. 1998). The best fitting amplitude A 0 b s is measured 
as explained above. To calculate A P g the PS of the dis¬ 
tribution is needed. In this study tracers are contained 
in simulation cubes with periodic boundaries, and the 
PS is obtained by Fourier transforming the distribution. 
Extrapolations of the PS beyond the observed fc-window 
necessary to calculate the integral for A P g introduce only 
a small error in A P g when d/y/2 < Rg < r max /10 where 
Rg is the Gaussian smoothing length, d is the mean 
tracer separation, and ?’ max is the sample size. 

3.3. Cluster/Void Abundance Parameters 

The third statistic is the numbers of clusters and voids. 
The statistic Ac is defined as 



A = 1 ) 3/2 

(2 tt) 2( 3 ’ 


G obs (v)dv/ / G bt {v)dv. 


( 8 ) 
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TABLE 1 

Characteristics of Cosmological N-Body Simulations 


Name 


U A 

h 

b 

N a 

N p 

L h 

Zi 

Nstep 

code 

LCDM1024 

0.27 

0.73 

0.71 

1.11 

2048 3 

2048 3 

1024 

17 

680 

PMTree 

LCDM5632 

0.27 

0.73 

0.71 

1.11 

2048 3 

2048 3 

5632 

17 

170 

PMTree 

LCDM410 

0.27 

0.73 

0.71 

1.11 

2048 3 

1024 3 409.6 

47 

470 

PM 

SCDM1024 

1.00 

0.00 

0.5 

1.5 

2048 3 

1024 3 

1024 

23 

230 

PM 

SCDM410 

1.00 

0.00 

0.5 

1.5 

2048 3 

1024 3 409.6 

47 

470 

PM 


“Size of mesh on which initial conditions are defined. 
6 (7i- 1 Mpc) 


where the integral is limited to 1.2 < v < 2.2 which is 
roughly centered at v = V3 which is the minimum in 
the best fit random phase curve where the number of 
isolated clusters is greatest. Ay is similarly defined over 
—2.2 < v < —1.2 roughly centered at v = — \/3, the 
minimum in the best fit random phase curve where the 
number of isolated voids is greatest. Here the best-fit 
Gaussian curve Gat is obtained as explained above using 
the middle part of the genus curve. The statistics Ac 
and Ay quantify the asymmetry of the genus curve and 
the multiplicity of clusters and voids. 

4. SIMULATIONS 

In order to study the behaviour of the genus statistic 
accurately we have made a set of large N-body simula¬ 
tions using a PM and PM+Tree codes. Our PM code 
uses the triangular-shaped cloud method for mass as¬ 
signment and force interpolation at particle positions. 
The force calculation at mesh points is done by the four- 
point finite difference algorithm (Park 1990, 1997). The 
PM+Tree code GOTPM, which is a merged version of 
our PM code with Dubinski’s tree code, is a parallel code 
which adds the tree force to the PM force at separations 
shorter than 4 times the pixel size of the mesh used in 
the PM part (Dubinski et al. 2003). For topology study 
of the matter distribution, the PM code is preferred since 
this code can make simulations with large dynamic range 
in mass at relatively low cost. The PM code is also pre¬ 
ferred when the peaks in the initial density field are used 
as biased tracers because high-resolution initial density 
fields are adopted in this code. Due to its high force res¬ 
olution power the PM+Tree code should be the choice 
when dark halos are targets for topology study. Because 
a single simulation cannot satisfy all of our needs we have 
made several simulations supplementing one another us¬ 
ing the PM and PM+Tree codes. Parameters of our sim¬ 
ulations are summarized in Table 1. 

Our results are mainly derived from PM+Tree simula¬ 
tions of the ACDM model with the WMAP parameters, 
namely, the density parameters fl m = 0.27, Da = 0.73, 
the Hubble constant ft = 0.71, <T m ,8 = 0.9, and the mat¬ 
ter is dominated by the cold dark matter (CDM). Here 
cr ni! 8 is the RMS mass fluctuation in a sphere of radius 
of 8 ft -1 Mpc. The physical sizes of the two simula¬ 
tion cubes are L = 1024 and 5632ft -1 Mpcs. We call 
them LCDM1024 and LCDM5632. We use the trans¬ 
fer function of the power spectrum given by Bardeen et 
al. (1986) for LCDM1024. But we have used Eisenstein 
& Hu (1998)’s power spectrum with Dbaryon = 0.0463 


for LCDM5632 to include baryonic effects. The num¬ 
ber of CDM particles evolved in both simulations is 
2048 3 = 8.59 billion, among the largest N-Body simu¬ 
lations ever made. A 2048 3 mesh is used for the PM 
part, and the force softening length in the Tree part is 
one tenth of the mean particle separation. Therefore, the 
dynamic range in length is 2 x 10 4 , from 0.05 to 1024 ft -1 
Mpc while the dynamic range in mass is 8.6 x 10 9 , from 
9.4 x 10 9 to 8.0 x 10 19 ft -1 M 0 in the case of LCDM1024, 
for example. Dark halos are identified in the LCDM1024 
simulation made by the PM+Tree code. Halo centers 
are found as local density peaks defined on a fine mesh. 
The CDM particles belonging to virialized dark halos are 
searched by checking their binding energy to the local 
halo centers and tidal radii of subhalos with respect to 
adjacent bigger halos (Kim & Park 2005). These iden¬ 
tification criteria overcome many of the problems and 
artificial results obtained from existing group-finding al¬ 
gorithms, and let us obtain a catalog of self-bounci phys¬ 
ical dark halos. Only the dark halos with more than 53 
member particles or heavier than 5 x 10 11 ft -1 M 0 are 
adopted. 

Even though our group finding algorithm allows us to 
identify many subhalos contained in big halos, those ha¬ 
los with mass much higher than 1O 13 M 0 cannot be con¬ 
sidered to be a single galaxy. On the other hand, every 
halo identified in the N-body simulation might not con¬ 
tain a galaxy. Therefore, even though dark halos might 
be better tracers of the observed galaxies than matter, 
they are still biased tracers. We address this problem 
by introducing the Halo Occupation Distribution (HOD) 
prescription to locate ‘galaxies’ within our dark halos. 
We adopt the following recipe of Zehavi et al. (2004). A 
central galaxy is assigned to a halo if its mass M exceeds 
Mmin, and the central galaxy can have satellites. The 
mean number of satellite galaxies is given by a power- 
law (-/Vgat) = (M/Mi)“, and they have a Poisson distri¬ 
bution. We use, as a toy model, the HOD parameters of 
logM m in = 11.76, logMi = 13.15, and a = 1.13 given by 
Table 3 of Zehavi et al. corresponding to galaxies with 
an absolute r magnitude cut of M r = —19.5 in the Sloan 
Digital Sky Survey. These ‘galaxy’ tracers are used in 
subsection 5.3. 

As an another galaxy formation model we adopt the 
biased galaxy formation scenario where galaxies form at 
peaks in the initial density field (Bardeen et al. 1986). 
We use PM simulations of 1024 3 CDM particles in a 
2048 3 mesh spanning 409.6ft. -1 Mpc along a side. In the 
ACDM model a set of peak bias parameters, v th = 0.05 
and R s = 0.763ft -1 Mpc, gives us peak particles which 
cluster like bright galaxies at the present epoch, where 
zz t h is the threshold density level for the peaks, and R s 
is the Gaussian smoothing length corresponding to the 
galaxy scale. In the case of the SCDM model, the peak 
parameters adopted are vta = 1-7 and R s = 0.585ft -1 
Mpc for cr mass> 8 = 0.6. These peak particles have been 
shown to follow the centers of collapsed dark halos, but 
are too clustered in high density regions since no merging 
process is included (Park & Gott 1991). 

5. RESULTS 

We have measured the genus statistics in the four di¬ 
mensional space of smoothing scale, redshift, tracers, and 
cosmogony. We use Gaussian smoothing and the defini- 



4 


Park et al. 



I_i_i_i_i_i_ i i i I 

1 10 
R g (h -1 Mpc) 


Fig. 1.— The genus-related statistics for dark matter measured 
from the LCDM1024 simulation. Dots with different sizes represent 
the statistics at redshifts z = 0,1, 2, 5, and 8. The largest dots 
connected by lines are for z = 0, and the smallest are for z = 8. 
In the bottom panel filled circles are Aq' s, and open circles are 
Ay's which are slightly shifted to avoid confusion. Stars are the 
statistics measured in redshift space. In the bottom panel bold 
stars are Ac's, and thin stars are Ay’ s. 


tion of the smoothing length Rq is standard (y/2 times 
larger than the ‘e-folding’ smoothing adopted by Gott et 
al. (1989)). We always maintain the condition that the 
Gaussian smoothing length Rq must be greater than or 
equal to the mean particle/halo separations. Topology 
is explored at 14 smoothing scales from Rq = 1.5ft -1 
Mpc to 30 ft -1 Mpc in the case of the matter distri¬ 
bution. Our large and high resolution simulations en¬ 
able us to obtain very accurate genus values. The num¬ 
ber of resolution elements in the simulation box is over 
2500 (= L 3 /(2tt) 3 / 2 Rq) at the largest smoothing scale 
Rg = 30ft -1 Mpc, and exceeds 2 x 10 7 at Rq = 1.5ft -1 
Mpc. The genus is also calculated at 7 epochs with 
redshifts 2 from 0 to 8, but we show results only at 
z = 0,1, 2, 5, and 8 to avoid over-crowding of points in 
figures. 

5.1. A CDM Dark Matter 

The mean particle separation of the LCDM1024 sim¬ 
ulation is 0.5 ft -1 Mpc, and we explore the topology of 
the dark matter distribution from Rq = 1.5ft -1 to 30 ft -1 
Mpc. The size of density array is 2048 3 , so the smooth¬ 
ing length is at least three times larger than the pixel 
size (thus avoiding discreteness effects). The results for 
the genus-related statistics are shown in Figure 1 as filled 
and open circles. Dots with different sizes represents the 
statistics at different redshifts from z = 0 (largest circles 
connected by solid lines), 1, 2, 5 to 8 (smallest). The dark 
matter in the ACDM universe clearly shows an intriguing 
development of small-scale topology in the cases of Av 
and Ay statistics in a direction opposite to that at large 
scales. The shift parameter Av rapidly goes to positive 
values at Rq < 9ft -1 Mpc. But at scales Rg > 9ft -1 
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Fig. 2.— Same as Figure 1, but statistics are calculated from 
the LCDM5632 simulation. Note that the data at z = 8 and 5 are 
almost overlapping. Only measurements in real space are plotted. 


Mpc it slowly moves in the negative direction and con¬ 
tinues to be negative even at 30 ft -1 Mpc scale. All these 
genus curves have a spongelike topology (G > 0) at the 
median density contour, so we would not call the negative 
and positive shifts of the central part of the genus curve 
as in the direction of ’meat-ball’ and ’bubble’ topologies, 
respectively, since their correspondence does not seem to 
be so simple. The transition scale of Av from positive 
to negative shift agrees with the analytic prediction of 
subsection 5.6. 

In Figure 2 the genus-related statistics calculated from 
the LCDM5632 simulation are plotted to show their be¬ 
haviour at larger scales. Note that the power spectrum of 
LCDM5632 is slightly different from that of LCDM1024, 
thus making the transition scale of Av occur at different 
Rg- It can be seen that non-Gaussianity is visible up 
to the scale ~ 100ft -1 Mpc. But using the analytic es¬ 
timation of the non-Gaussianity effects of gravitational 
evolution in this weakly non-linear scales given in sub¬ 
section 5.6, one could draw conclusions on the Gaussian- 
ity of the initial density field from data with smoothing 
lengths much smaller than 100ft -1 Mpc if the matter dis¬ 
tribution can be used. 

The parameter Ay for voids (open circles in the bottom 
panels of Figure 1 and 2) shows a behaviour similar to the 
shift parameter. At scales below 10 ft -1 Mpc voids are 
broken into many pieces and the number of voids exceeds 
the Gaussian value predicted by the middle part of the 
genus curve (see Figure 1). But at larger scales voids 
percolate and the number of voids is slightly less than 
the Gaussian value (see Figure 2). On the other hand, 
the parameter Ac for clusters (filled circles in the bottom 
panel of Figure 1 and 2) behaves differently. It tends to 
increases at all scales at the present epoch. The speed 
with which new high density clumps stick up above the 
threshold and appear in the sample is faster than that of 
merging and disappearance of existing clumps. But at 
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very small scales (Re < 2 ft -1 Mpc) the merging speed 
wins over the sticking-up speed at high redshifts (z > 5) 
and subsequently the situation reverses. At large scales, 
the sticking-up of clusters and merging of voids are the 
dominant processes (Figure 2). It should be noted that 
gravitational evolution only makes the both parameters 
Ay and Ac increase over 1 at scales below 10 ft -1 Mpc at 
the present epoch. Therefore, the observed value of the 
void abundance parameter Ay which is less than 1 (Park 
et al. 2005) must be due to strong biasing effects, which 
tend to make small voids merge with one another. This 
can be achieved if galaxies do not form in weak walls and 
filaments of large scale structure dividing small voids. 

The amplitude drop parameter shown in the top panel 
of Figure 1 evolves very rapidly at scales below about 10 
ft -1 Mpc as the gravitational evolution proceeds. And its 
scale dependence is again determined by the characteris¬ 
tic shape of the power spectrum of the ACDM model. 

Also shown in the figures are the genus-related statis¬ 
tics measured in redsliift space at z = 0 (stars in all 
figures). We adopt the far-field approximation to the 
redslrift distortion of the density field. Displacement of 
particles is thus assumed to occur only along the direc¬ 
tion of a coordinate axis. We average three genus curves 
measured from redshift space density fields each of which 
is perturbed in the direction of each coordinate axis using 
the periodic boundary condition. Genus and the power 
spectrum of density field in redshift space are measured 
to calculate the genus-related statistics. The amplitude 
drop parameter Ra increases in redshift space because 
the amplitude of the random-phase genus curve (Aps in 
Equation (7)) drops more than that of non-linear genus 
curve (A 0 bs) due to the randomization effects of small- 
scale peculiar velocity. Effects of redshift space distor¬ 
tions on the genus curve have been previously studied 
by Melott, Weinberg, & Gott (1988). They have re¬ 
ported that the effects are negligible when a smoothing 
larger than the correlation length is applied to the den¬ 
sity field. And the linear regime calculation of Matsubara 
(1996) has supported their claim. Subsequently, redshift 
space distortion effects have been generally ignored in the 
genus analyses. This is qualitatively true for the matter 
field of the ACDM model where the correlation length 
is about 5.5 ft -1 Mpc a,t z = 0 and smoothing lengths 
larger than this are studied. However, as can be seen in 
our figures, there are non-negligible systematic effects of 
redshift space distortion on the genus related statistics 
at the small smoothing scales we are interested in. An 
example is the Ac parameter (thick stars in the bottom 
panel of Figure 1). Being a measure of the topology of 
high density regions, Ac is most sensitive to the peculiar 
velocity field. 

5.2. ACDM Biased Peak Particles 

We use the peaks in the initial Gaussian density field 
as the biased ‘galaxy’ particles. They are defined by 
the threshold level vth and the galaxy scale R s (Park & 
Gott 1991), which are chosen so that the resulting bi¬ 
ased particles approximate the distribution of galaxies 
at the present epoch. When we require them to have 
°pk,8 = 1 and to have a correlation function close to the 
observed galaxy CF at z = 0, we find a set, Pth = 0.05 
and R s = 0.763 ft -1 Mpc in our ACDM universe with 
Cm,8 = 0.9. The number density of the peaks is much 
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Fig. 3.— Genus-related statistics for the biased peaks calculated 
from the ACDM410 simulation. Notations are the same as Figure 
1. 


lower than that of CDM particles, and the genus is mea¬ 
sured only at scales Rq > 5ft -1 Mpc. The simula¬ 
tion used for identifying the biased peak particles is the 
ACDM model evolved by the PM code in a 2048 3 mesh 
spanning 409.6ft. -1 Mpc along a side. Its small pixel 
size (0.2ft -1 Mpc) allows us to look for density peaks 
at the galaxy scale R s . The genus-related statistics of 
the biased peak particles are shown in Figure 3. Even 
though our biased density peaks collapse at a specific 
epoch and can be considered as galaxies after then, we 
plot the statistics of the biased particles at all epochs 
after z = 8 for theoretical interest. 

They behave essentially in the same way as those for 
cold dark matter, but there are some changes. Similarly 
to before, the peak particles develop a negative shift Av 
at scales larger than 6ft -1 Mpc, and a positive shift at 
smaller scales. And the cluster abundance parameter 
Ac increases at all scales above 5ft -1 Mpc. However, 
the void abundance parameter Ay (open circles in Fig¬ 
ure 3) shows more changes. Namely, at scales smaller 
than about 7 ft -1 Mpc the biasing prescription makes 
small voids percolate and Ay become less than 1 at high 
redshifts. As noted above, gravitational evolution only 
makes Ay increase above 1 at scales less than 10 ft -1 Mpc, 
contrary to observations. Here we have a specific exam¬ 
ple of a galaxy formation mechanism which allows exist¬ 
ing voids to percolate strongly enough to overcome the 
nonlinear gravitational effect of opening of new voids and 
makes Ay less than 1. 

The stars in Figure 3 are the statistics measured in 
redshift space at 2 = 0. The effect significantly decreases 
Ac (bold stars in the bottom panel of Figure 3) at all 
scales explored as dense clusters are spread along fingers- 
of-gocl. It can be noted that the void abundance param¬ 
eter is not much affected by the redshift space distortion 
because empty regions are much more stable to the dis¬ 
tortion compared to dense clumps. This observation is 
also true for other tracers used in this work. Given its 
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insensitivity to redshift space distortion, Ay is a crucial 
statistic telling the story of galaxy formation in a less 
distorted way. 

5.3. A CDM Dark Halos 

The biased particles studied above are the density 
peaks in the Gaussian initial density field. This struc¬ 
ture formation scheme assumes that galaxies form at the 
maxima in the primordial density fluctuations. It can be 
shown that these actually correspond well to collapsed 
dark halos at the final epoch of simulations (Park & 
Gott 1991). However, it would be good to check the idea 
and to re-compute the same statistics by using physically 
bound dark halos identified in high resolution simula¬ 
tions. Another important difference between the biased 
peak particles and dark halos is that formation of tracers 
is naturally incorporated in the genus analysis in the case 
of dark halos. Since we need high resolution for halo for¬ 
mation and identification and large volume for high sta¬ 
tistical significance, the LCDM1024 simulation is used. 
The minimum halo mass identified is 5 x 1O 11 M 0 (53 
particles). With the dark halo samples obtained from 
this simulation, topology of the dark halo distribution is 
studied at scales from 5 to 30 ft -1 Mpc at z = 0. At 
higher redshifts the minimum smoothing length used be¬ 
comes larger as the dark halos collapsed by then become 
rarer. 

The genus-related statistics measured for the dark ha¬ 
los in the ACDM model are shown in Figure 4. Topology 
of ACDM dark halos is very different from those of under¬ 
lying dark matter or biased peak particles. The shift pa¬ 
rameter Au starts from large negative values at high red- 
shifts when the only densest regions can collapse. Con¬ 
sistent with this behaviour the Ac parameter is greater 
than 1, and Ay is lesser than 1 at high redshifts. At 
lower redshifts, the topology of the dark halo distribution 
evolves strongly. Ay becomes positive and Ac becomes 
less than 1. This is opposite to what we have seen in 
the case of CDM matter and biased peak particle distri¬ 
butions. Ay still stays below 1 but has increased from 
its high redshift value, which is consistent with the case 
of biased peak particles. Ac of dark halos rapidly de¬ 
creases at low redshift while that of peak particles keeps 
increasing. This is because the dark subhalos merge to 
form one huge halo at centers of clusters of halos, and 
thus decreasing the number of subhalos or reducing the 
richness of clusters as non-linear evolution proceeds at 
small scales. This merging process within clusters does 
not exist for peak particles because peak particles never 
lose their identity and the clusters of peak particles get 
richer as they become more compact. The reality might 
be somewhere between these two extremes. Comparison 
of Figure 4 with previous ones demonstrates sensitivity of 
the topology to the galaxy formation mechanism. Since 
galaxies show different topology depending on their in¬ 
ternal physical properties, topology analysis can be used 
to discriminate among different galaxy formation mech¬ 
anisms for different species of galaxies. 

Hikage, Taruya & Suto (2003) have calculated the bi¬ 
asing and genus of dark matter halos as a function of 
halo mass and smoothing scale by using the Hubble Vol¬ 
ume Simulation. They have found that the halo bi¬ 
asing effects on the genus are comparable to the non- 
Gaussianity due to the non-linear gravitational evolu- 



Fig. 4.— Genus-related statistics for the dark halos calculated 
from the LCDM1024 simulation. Notations are the same as Figure 
1. Since the number density of dark halos at z = 8 and 5 is too 
low, we show the results at z = 0,1, 2, and 3.5. And at z = 3.5 
(smallest points) genus is measured at Rq > 7 h ~ 1 Mpc. 


tion. This is consistent with our results that changes 
of the genus-related statistics due to gravitational evo¬ 
lution at Re = 30ft -1 Mpc scale is of the order of that 
due to tracer difference. Hikage et al. (2003) have stud¬ 
ied topology of halo distribution at Re > 30ft -1 Mpc. 
However, at Re = 30ft -1 Mpc their results are thought 
to be affected by discreteness effects because the mean 
separation of their halo subsets is only about 44 ft -1 
Mpc. In topology analysis of a set of points, the smooth¬ 
ing length should be larger than the mean tracer sepa¬ 
ration to avoid systematic discreteness effects, or some 
correction should be made to remove the effects (Park 
et al. 2005). The large pixel size of 23ft -1 Mpc used 
in their analysis may have also caused some discreteness 
effects at that smoothing scale. In general one wants the 
smoothing length Re to be at least 1.77 times the pixel 
size to avoid discreteness effects (eg. Hamilton, Gott, & 
Weinberg 1986). To identify halos they have used the 
friends-of-friends algorithm, which is known to miss sub¬ 
halos within big halos and often produces unphysical ha¬ 
los connected by narrow chains of particles. Therefore, it 
is not easy to directly compare their results with those of 
the current work in detail. But our results qualitatively 
agree with their conclusion that non-Gaussianity induced 
by halo biasing is comparable to that by the non-linear 
gravitational evolution. 

5.4. HOD Galaxies 

A halo identified in the N-body simulation might not 
contain a galaxy, or can host more than one galaxies. 
Therefore, even though dark halos might be a better 
tracer of the observed galaxies than matter, they may 
not be representing galaxy distribution well enough to 
allow quantitative comparisons. We address this prob¬ 
lem by introducing the Halo Occupation Distribution 
(HOD) prescription to locate ‘galaxies’ within our dark 
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Fig. 5. — Genus-related statistics for the ‘galaxies’ derived from 
the dark halos of the LCDM1024 simulation. Notations are the 
same as Figure 1. Only the results at z — 0 are shown. 
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Fig. 6.— Genus-related statistics for dark matter measured from 
the SCDM1024 simulation. Notations are the same as Figure 1. 


halos. We adopt the following recipe of Zehavi et al. 
(2004). A central galaxy is assigned to a halo if its mass 
M exceeds M m j n , and the central galaxy can have satel¬ 
lites. The mean number of satellite galaxies is given by a 
power-law (N sat ) = and they have a Poisson 

distribution. The satellites are randomly located at dis¬ 
tances of 0.3ft. _1 Mpc from the central galaxy. This does 
not change our results because the minimum smooth¬ 
ing scale adopted for analysis of the HOD galaxies is 
5fo _1 Mpc. We use, as a toy model, the HOD parameters 
of logM m i n = 11.76, logM-i = 13.15, and a = 1.13 given 
by Table 3 of Zehavi et al. corresponding to galaxies 
with an absolute r magnitude cut of M r = —19.5 in the 
Sloan Digital Sky Survey. 

In Figure 5 we show the genus-related statistics cal¬ 
culated from distribution of ‘galaxies’ assigned within 
dark halos of the LCDM1024 simulation according to the 
HOD prescription. The correlation length of these model 
galaxies is 5.6/i -1 Mpc, similar to that of bright optical 
galaxies. Only the results at z = 0 are shown. Compar¬ 
ing the statistics with those of dark halos at 2 = 0, we 
find Av becomes less positive, Ac is increased, and Ay 
is decreased. These changes the topology of our HOD 
‘galaxies’ to be more compatible with that observed for 
SDSS galaxies (Park et al. 2005). Topology depends on 
the type of observed galaxies, but one phenomenon con¬ 
sistently found in observations is that the Ay parameter 
is less than 1. And This behaviour is seen in the distri¬ 
bution of our HOD ‘galaxies. However, the observed Ac 
parameter is typically greater than Ay, which does not 
quite agree with the HOD ‘galaxies’. The HOD method 
seems to be promising in explaining not only the one- 
point and two-point functions of galaxy distribution, but 
also functions of high-order moments like genus. How¬ 
ever, more detailed comparisons have to be made to con¬ 
firm it. 

Benson et al. (2001) has used a semi-analytic model 


of galaxy formation to obtain the distribution of model 
galaxies, and to measure the genus curves at smoothing 
lengths of 4.24 and 5.66 ft, -1 Mpc at redshifts z = 0 and 
1. But they have not been able to detect the difference 
between the genus curves for the dark matter and the 
model galaxies due to sparseness and smallness of their 
simulation data. 

5.5. SCDM Case 

As a fiducial model we consider the SCDM model to 
study the dependence of topology on cosmology. In this 
model = 1 and h = 0.5, and we adopt a bias factor 
b = 1.5 and Bardeen et al. (1996)’s fitting formula for the 
power spectrum. Figure 6 shows that the matter distri¬ 
bution of the SCDM universe has qualitatively the same 
topology as that of the ACDM universe. Only the char¬ 
acteristic scales are moved down in the cases of Av, and 
Ay . The switching from negative to positive shift occurs 
at ~ 4 h~ x Mpc, and the switching from sub-Gaussian to 
super-Gaussian amplitude for the void part of the genus 
curve occurs at a similar scale. This change in the char¬ 
acteristic scales is the direct consequence of the change in 
the shape of power spectrum. Figure 7 shows the statis¬ 
tics for the biased peak particles found in the SCDM sim¬ 
ulation. For Av a negative shift prevails at all scales, and 
the redshift space distortion effects make little change in 
Av. Gravitational evolution makes the Ac parameter 
increase above 1 as in the ACDM case while the red- 
shift space distortion effects make it decrease. The void 
abundance parameter Ay is significantly less 1 due to 
the biasing effects at small scales, and is hardly affected 
by the redshift space distortion as in the ACDM case. 

5.6. Comparison with Perturbation Theory 

Matsubara (1994) has obtained a formula for the genus 
curve modified due to gravitational evolution in the 
weakly non-linear regime by perturbatively expanding 
the statistic in ctq = (<5p/p) rms , the RMS fluctuation of 
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Fig. 7.— Genus-related statistics for the biased peaks measured 
from in the SCDM410 simulation. Notations are the same as Figure 
1. 


an over-density field smoothed over a given scale. The 
non-linear correction in the genus curve to the first order 
in do is an odd function of v and, therefore, causes a shift 
and an asymmetry between the high and low density re¬ 
gions. It is hoped that the gravitational evolution effects 
on genus in the weakly non-linear regime can be mod¬ 
eled by this analytic theory. When the volume threshold 
level is used, the three-dimensional genus expanded to 
the first order in cr 0 is given by (Matsubara 2003) 


G(y f ) = 







+ (S (2) _ SW)Vf](To), 

(9) 

where the variance parameters are 
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and the window function is a Gaussian W(kR) = 
exp(— k 2 R 2 /2). SM are the skewness parameters defined 
as 
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where k± 2 = |ki + k 2 | and B is the bispectrum. The 
skewness parameters can be calculated from (Matsubara 
2003) 
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Here P is the linear power spectrum and I v are the mod¬ 
ified Bessel functions. 

We have performed the above integrals to get the skew¬ 
ness parameters for our ACDM and SCDM models as a 
function of smoothing length Rg- The curves in Figure 
8 show the genus-related statistics predicted by equa¬ 
tion (9) for the ACDM and SCDM models at the present 
epoch z = 0. In the top panel the RMS density fluctu¬ 
ation do and the skewness parameters are shown. The 
differences between the skewness parameters in equation 
(9) are very small as can be seen here, and make the 
genus as a function of Vf rather insensitive to gravita¬ 
tional evolution. According to the perturbation theory 
the amplitude of the genus curve remains the same in 
the weakly non-linear regime. This is certainly not true 
because the density field starts to build phase correla¬ 
tions and the genus amplitude starts to decrease as soon 
as the gravitational evolution occurs. The behaviours of 
the At/, Ac, and Ay parameters predicted by equation 
(9) are remarkably similar to those seen in Figure 1 and 
6. This is an important confirmation of the Matsubara’s 
theory. The one major failure of the analytic prediction is 
for the Ac parameter at small non-linear scales. At high 
redsliifts Ac decreases below 1 as high density clumps 
merge together as predicted by the analytic theory. But 
at lower redshifts Ac increases to become greater than 
1 as the merging rate slows down and the speed new 
clumps appear becomes higher. 

Matsubara (1996) has derived the genus in redshift 
space in the linear regime. The genus in redshift space 
predicted by the linear theory is 


G {s) H = ^VC(1 - C)G {r) (v), 
where G(y) is the genus curve in real space, 


(19) 
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Fig. 8. — Genus-related statistics at z = 0 according to the 
Matsubara’s perturbation theory. In the second and third panels 
thick curves are for scales with ctq < 0.25, medium thickness curves 
for 0.25 < cto < 0.5, and thin curves for 0.5 < ao < 1.0. In the 
top panel the RMS deviation of density fluctuation and skewness 
parameters are plotted. 


b is the bias factor, and 
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where D is the linear growth factor, and a is the ex¬ 
pansion parameter. The formula says the redshift space 
distortion does not alter the shape of the genus curve 
and only affects its amplitude. For matter in our ACDM 
model b = 0.9, / = 0.468 and C = 0.414, and the ampli¬ 
tude drops by a factor of 0.980 while in the SCDM model 
which has b = 1.5, / = 1, and C = 0.433 the factor is 
0.970. In the linear regime the genus amplitude calcu¬ 
lated from redshift space power spectrum is equal to that 
from real space power spectrum since the shape of power 
spectrum remains the same. This makes the denomina¬ 
tor of the amplitude drop parameter Ra unchanged in 
redshift space, and makes Ra itself decrease by above 
factors in the ACDM and SCDM models. This can be 
observed in Figures 1 and 6 at scales Rq > 20/i _1 Mpc. 

Matsubara & Suto (1996) have used N-body simula¬ 
tions to examine the redshift space distortion effects on 
genus. They have reported that the amplitude of the 
genus curve is suppressed more than that expected by 
linear theory, but that its shape in redshift space still 
remains the same down to weakly non-linear scales. We 
have also shown that the redshift space distortion effects 
on the genus curve are small in the linear and quasi-linear 
scales. But it should be pointed out that the high thresh¬ 
old part of the genus curve (namely the Ac parameter) 
is more affected. 

Comparisons of numerically calculated genus with the¬ 
oretical predictions are useful. But so far theories have 
mainly concerned themselves with the development of 
non-Gaussianity in the local density field rather than 
changes in topology. Our genus meta-statistics are di¬ 
rectly addressing questions about the gravitational evo¬ 


lution of topology and the effects of biasing and redshift 
space distortion on topology. 

6. CONCLUSIONS 

We have studied dependence of the genus statistic 
on cosmology, tracer, redshift, and Gaussian smoothing 
scale. We summarize our findings. 

1. The shift parameter A^ of the genus curve shows a 
strong scale dependence whose characteristics are deter¬ 
mined by the shape of the initial power spectrum. For the 
matter distribution in the ACDM model the negative-to- 
positive transition occurs at ~ 9h -1 Mpc while it is at 
about 4 h -1 Mpc for SCDM. Below these scales grav¬ 
itation evolution causes strong positive shifts starting 
from high redshifts. At larger scales the gravitational 
evolution produces weak negative shifts. This trend en¬ 
hances as the evolution proceeds. When the topology 
tracers are biased peak particles instead of CDM parti¬ 
cles, evolution of the genus curve is qualitatively simi¬ 
lar. But the positive-to-negative transition scale is dif¬ 
ferent, and at small scales A v keeps slowly decreasing, 
which is opposite to the case of CDM matter particles. 
We have demonstrated that galaxies can have a topol¬ 
ogy very different from that of matter on small scales. 
The Av parameter measured from the dark halo and 
HOD ‘galaxy’ distributions is positive at scales between 
5 and 30 /i _1 Mpc, and becomes more positive due to the 
combined effects of formation of new objects and grav¬ 
itational clustering. This trend is opposite to that seen 
in distributions of matter and biased peak particles at 
scales larger than 9h -1 Mpc. It remains to be seen which 
value of the shift parameter the distribution of real galax¬ 
ies has. The redshift space distortion systematically de¬ 
creases A v in the ACDM universe, and the effects are 
significant at scales < 5/i -1 Mpc. We can in principle 
correct for these redshift space distortion effects on Av 
by using the observationally measured line-of-sight ve¬ 
locity dispersion to correct the shift parameter at both 
large and small scales for the study of primordial density 
field and galaxy formation mechanisms. 

2. The amplitude drop parameter is more difficult to 
measure than the genus curve alone because it requires 
additional accurate measurement of power spectrum for 
each tracer of topology. Direct consequence of genus am¬ 
plitude drop is the phase correlation of density held, and 
the phase correlation depends on the shape of the ini¬ 
tial density power spectrum and the biasing. We have 
found that the dark halos and the HOD ‘galaxies’ show 
a very small genus amplitude drop with 0.9 < Ra < 1 
down to the scale 5/) _1 Mpc with little scale dependence. 
And it hardly changes due to gravitational evolution or 
peculiar motions. This might be incompatible with the 
observationally measured value of about 0.6 at 5/i -1 Mpc 
for the IRAS Point Source Catalogue Redshift Survey 
(Canavezes et al. 1988) or for the Center for Astrophysics 
Survey (Vogeley et al. 1994). Satisfying the constraint 
from the genus amplitude drop can be a challenge for 
HOD modeling of galaxy formation. 

3. The void abundance parameter Ay of the dark 
matter distribution strongly increases above the Gaus¬ 
sian value below 10 h _1 Mpc, but gently drops below 1 
at larger scales as the gravitational evolution proceeds. 
This is true for both ACDM and SCDM models even 
though their transition scales are different. While grav- 
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itational evolution always raises Ay at small scales, bi¬ 
asing can lower it below 1 as can be seen for the biased 
peak particles, dark halos, and HOD ‘galaxies’. Even the 
direction of evolution of Ay (and Ac) can be changed by 
biasing. Since Ay is observed to be less than 1 for the 
SDSS galaxies (Park et al. 2005), this fact is strong ev¬ 
idence for the existence of biasing in the distribution of 
galaxies with respect to matter. Fortunately, Ay at small 
scales is hardly affected by the redshift distortion effects, 
and can give faithful constraints on galaxy formation in 
under-dense regions. 

We have considered only three mechanisms for galaxy 
formation. The peak biasing scheme is able to explain the 
stronger clustering amplitude for rarer objects (Kaiser 
1984). But we have considered only one version of the 
peak biasing scheme where the threshold for galaxies is 
a step function, and have not taken into account merg¬ 
ing (cf. Narayanan, Berlind, & Weinberg 2000 for var¬ 
ious biasing models). The first assumption will affect 
the voids, and the second the clusters. The arbitrari¬ 
ness of the threshold function can be removed by using 
the dark halos which have collapsed and are self-bound 
physical objects. But the dark halos also cannot fully 
trace galaxies because the mass function of dark halos 
continuously extends up to supercluster scales and does 
not show a natural cutoff for galaxies. Some of them 
with small masses may have lost baryons in high den¬ 
sity environments, and have not been able to form stars. 
The HOD ‘galaxies’ seem to approximate the distribu¬ 
tion of real galaxies better than the dark halos because 
the Ay parameter is closer to observations for the HOD 
‘galaxies’ even though the cluster abundance parameter 
Ac may be still too small (cf. Park et al. 2005). 

While galaxy formation and evolution in a high den¬ 
sity environment are complicated due to the merging of 
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structures and the high ram pressure of the intergalac- 
tic medium (Gunn & Gott 1972), it may be relatively 
easier to model galaxy formation in low density regions. 
For example, biased peak particles, dark halos, and HOD 
‘galaxies’ all produce Ay smaller than 1 consistent with 
observations. But their Ac are wildly different. There¬ 
fore, even though identification of the biased peaks and 
dark halos with galaxies is not a good approximation in 
high density regions, they seem to trace galaxies better 
in under-dense regions for which there is a sensitive mea¬ 
sure like the void abundance parameter Ay. It should 
be emphasized that the cosmic voids have not been stud¬ 
ied as extensively as clusters and superclusters, and it is 
only in recent years that voids have begun to be actively 
studied (Hoyle & Vogeley 2004; Sheth & van de Weygaert 
2004; Rojas et al. 2004; Colberg et al. 2004). After all, 
voids are not the places where there is nothing, but places 
where history of the universe is better kept. With the ad¬ 
vent of new large observational data sets like the SDSS 
sample we can hope better to understand the formation 
mechanism of each species of galaxies through the topol¬ 
ogy analysis. We are currently studying the topology 
of galaxies divided into subgroups with different internal 
physical properties. 
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